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Abstract 

Plasmodium of Physarum polycephalum is a large cell capable of solving graph- 
theoretic, optimization and computational geometry problems due to its unique 
foraging behavior. Also the Plasmodium is unique biological substrate that mimics 
universal storage modification machines, namely the Kolmogorov-Uspensky ma- 
chine. In the Plasmodium implementation of the storage modification machine data 
are represented by sources of nutrients and memory structure by protoplasmic tubes 
connecting the sources. In laboratory experiments and simulation we demonstrate 
how the plasmodium-based storage modification machine can be programmed. We 
show execution of the following operations with active zone (where computation oc- 
curs): merge two active zones, multiple active zone, translate active zone from one 
data site to another, direct active zone. Results of the paper bear two-fold value: 
they provide a basis for programming unconventional devices based on biological 
substrates and also shed light on behavioral patterns of the Plasmodium. 

Keywords: Physarum polycephalum^ Kolmogorov-Uspensky machine, pattern forma- 
tion, morphogenesis, graph theory 



1 Introduction 



Physarum polycephalurrp~\ commonly known as a true or multi-headed slime 
mold, is — at one stage of its complicated life cycle - a single cell with many 
diploid nuclei. This syncytial mass of protoplasm, called Plasmodium, looks 
like amorphous yellowish mass. The Plasmodium behaves and moves as a gi- 
ant amoeba. It feeds on bacteria, spores and other microbial creatures. When 
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foraging for its food the plasmodium propagates towards sources of food par- 
ticles, surrounds them, secretes enzymes and digests the food. TypicaUy the 
Plasmodium forms a congregation of protoplasm in a food source it occupies. 
When several sources of nutrients are scattered in the Plasmodium's range, 
the Plasmodium forms a network of protoplasmic tubes connecting the masses 
of protoplasm at the food sources. When we think of food sources as nodes 
and protoplasmic tubes as edges, we say the plasmodium develops a planar 
graph. 

Nakagaki et al [16,17,18] showed that the topology of the Plasmodium's proto- 
plasmic network optimizes the Plasmodium's harvesting on distributed sources 
of nutrients and makes more efficient flow and transport of intra-cellular com- 
ponents. Therefore the plasmodium is considered as a parallel computing sub- 
strate complementary [5] to existing massive-parallel reaction-diffusion chemi- 
cal processors [2] . Experimental observations suggest that during development 
of its protoplasmic network the plasmodium undergoes transitions between 
various classes of proximity graphs [6]. It starts with spanning trees [3] and 
may complete its protoplasmic network development as a Delaunay triangu- 
lation [22]. 

Implementation of a general purpose computing machine is the most remark- 
able feature of the plasmodium of Physarum polycephalum. In [4] we experi- 
mentally demonstrated that the plasmodium can implement the Kolmogorov- 
Uspensky (KUM) machine [13,34], a mathematical machine in which the stor- 
age structure is an irregular graph. The KUM is a forerunner and direct 'an- 
cestor' of Knuth's linking automata [12], Tarjan's reference machine [29], and 
Schonhage's storage modification machines [23,24]. The storage modification 
machines are basic architectures for random access machines, which are the 
basic architecture of modern-day computers. The plasmodium-based imple- 
mentation of KUM [4] provides a first-ever biological prototype of a general 
purpose computer. 

The key component of the KUM is an active zone [13,34], which may be seen as 
a computational equivalent to the head in a Turing machine. Physical control 
of the active zone is of utmost importance because it determines functionality 
of the biological storage modification machine. 

In the present paper we show — in laboratory and computer experiments with 
Physarum polycephalum — how basic operations Add node, Add edge, 
Remove Edge are implemented in the Physarum machine. We also provide 
unique results on controlling movement of an active zone. The paper is struc- 
tured as follows. We provide a very brief introduction to Kolmogorov-Uspensky 
machine in Sect. 2. Section 3 describes a simple experimental setup for study 
of Plasmodium. A particle-based model of the plasmodium is presented in 
Sect. 4. In Sect. 5 we discuss results of laboratory experiments and computer 
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Fig. 1. An exemplar snapshot of Physarum machine. Protoplasmic tube connecting 
flakes a and b represents an edge e(a^b) of KUM. Active zone A(c) emerges in the 
node c. 

simulation on reconfiguration of basic protoplasmic graphs. Overview of the 
results and future experiments are outlined in Sect. 6. 



2 Physarum machine 



Kolmogorov-Uspensky Machine (KUM) is defined on a labeled undirected 
graph (storage structure) with bounded degrees of nodes and bounded num- 
ber of labels [13,34]. KUM executes the following operations on its storage 
structure: select an active node in the storage graph; specify local active zone, 
i.e. the node's neighborhood; modify the active zone by adding a new node 
with the pair of edges, connecting the new node with the active node; delete 
a node with a pair of incident edges; add/delete the edge between the nodes. 
A program for KUM specifies how to replace the neighborhood of an active 
node (i.e. occupied by an active zone) with a new neighborhood, depending 
on the labels of edges connected to the active node and the labels of the nodes 
in proximity of the active node [8] . 

Physarum machine is a biological implementation of KUM, where a node of 
the storage structure is represented by a source of nutrients (e.g. an oat fiake); 
an edge connecting two nodes is a protoplasmic tube linking two sources of 
nutrients corresponding to the nodes; and, an active zone is domain of space 
(which may include food sources) occupied by a propagating pseudopodium. 
An sample architecture of Physarum machine is illustrated in Fig. 1. 
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In Physarum machine the computation is implemented by active zone, or sev- 
eral active zones. To make the computation process programmable one needs 
to find ways of sensible and purposeful manipulation with the active zones. 
Several operations, manipulation procedures, are discussed in the present pa- 
per. 



3 Experimental method 

The Plasmodia of Physarum polycephalum were cultured on wet paper towels, 
fed with oat fiakes, and moistened regularly. We subcultured the Plasmodium 
every 5-7 days. Experiments on growing spanning trees were undertaken in 
standard Petri dishes, 9 cm in diameter. A substrate was a wet filter paper. 
We preferred the filter paper not 2% agar gel, because the paper offers less 
favoured conditions for the plasmodium growth, and thus less branching of 
the propagating pseudopodia is observed (see details in [7,28]). The Petri 
dishes with plasmodia were kept in darkness and only exposed to light during 
observation and recording of images. 

Data points, to be connected by protoplasmic graphs, were represented by oat 
fiakes. Photographs of the protoplasmic networks, developed by the Plasmod- 
ium, were made using the digital camera FinePix S6500. Colours of original 
images were enhanced by increasing contrast to 100 and decreasing gamma 
correction to 0.3. 



4 Particle-based model approximating the behaviour of Physarum 
Plasmodium 

To provide the computer model approximation of Physarum plasmodium ac- 
tive zone computation, the laboratory experiment configurations were used to 
map synthetic environments for the model organism. The model is based on the 
simple particle behaviours that were used in [10] to generate emergent trans- 
port networks. The model is computationally very simple, using distributed 
local sensory behaviours, but approximates some of the complex phenom- 
ena observed in Physarum: foraging for food stimuli, spatially distributed 
oscillations, oscillation phase shifting, shuttle streaming, amoebic movement, 
network formation, surface area minimisation, and network minimisation be- 
haviours. A full description of the model is in preparation, but an overview 
specifically in relation to plasmodium behaviour is given below. 

The model takes a multi-agent (particle) approach to generate emergent plas- 
modium behaviours. Movement and internal oscillations of the plasmodium is 
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Fig. 2. Problem data projected to coupled map storing diffusing food source and 
agent trails. 

based upon the collective behaviour of the agent population. The movement 
of agents corresponds to the flux of sol within the Plasmodium. Cohesion of 
the Plasmodium is an emergent surface tension property that arises due to the 
agent-agent interactions. Directional orientation and movement of the Plas- 
modium active front is generated by coupling the emergent mass behaviours 
with chemoattraction to local food source stimuli. 

A population of mobile agents, adopting simple stimulus-response behaviours 
is created. The population adopts autocrine chemotaxis behaviours — agents 
both secrete and sense approximations of chemical trails. The population is 
coupled to a two-dimensional discrete map representing the problem conflgu- 
ration — the 'data' map, Fig. 2. 

Food source stimuli whose positions are stored in the data map are projected 
to another coupled map (the 'trail' map) at every step of the scheduler. The 
strength of the projected food sources can be adjusted with a weighting param- 
eter. When the Plasmodium engulfs a food source the stimulus for diflFusion is 
reduced by the encapsulation. This is approximated in the model by damping 
the food stimulus projection to 1% of the food value if the site is occupied by 
an agent. The food stimuli are diffused by means of a simple 3x3 mean fllter 
kernel. The combined diffusion stimuli and agent trails stored in the trail map 
are degraded at every scheduler step by 0.01 to maintain a relatively steep dif- 
fusion gradient away from the food source. The diffusion gradient corresponds 
to the quality of the nutrient and substrate of the Plasmodium's environment 
(for example the different growth patterns seen in soaked fllter paper and agar 
substrates), and differences in the stimulus strength, stimulus area, affect both 
the steepness, and propagation distance of the diffusion gradient and affect the 
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Fig. 3. Morphology of agent particle, 
growth patterns of the synthetic Plasmodium. 

The agent population (size 2000 agents for all experiments by restructuring 
of the protoplasmic graph, in last case we employed 3200 agents) is initialised 
with random agent positions (on lattice 200 x 200 sites) and random agent 
orientations. To reproduce the 'wet' Physarum experiments agents were in- 
oculated at particular locations (for example food sources) and new agents 
are initialised at the front of Plasmodium growth according to local occu- 
pancy measures. These measures are based on a local window surrounding 
each agent {6yj)^ a minimum occupancy threshold (^min), and a maximum 
threshold (^max)- In these results = 3x3 window, ^min > and ^^ax < 4. 
The morphology of each agent is shown in Fig 3. Let /, 5, r be concentrations 
of chemo-attractants as measured by left (FL in Fig. 3), forward (F in Fig. 3, 
and right (FR in Fig. 3) sensors. Each agent can be characterised at time step 
t by its current position and the angle the agent rotates at this time step. 
At each step of the simulation the rotation angle a is calculated as follows: 



0, 

random(— a, a). 



if > and > 

if > and > 

if > and > 
otherwise 



Agent's position is updated as follows: 



R(a) o T{5,), if 7(R(a) o T(4)) = 

R(random(360)), otherwise 



where R(a) is rotation on angle a operation and T{Sx) is a translation on 
distance 6x operation, and for any Euclidean point y ^{y) = 1 if there is an 
agent residing at site y] ^{y) = otherwise. 
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The agents sense the concentration of stimuh in the trail map from three 
forward oriented sensors. At each scheduler step the agents orient themselves 
towards the strongest trail source by rotating left or right, depending on the 
source of the stimuh. 

For these simulations sensor width is one pixel, sensor arm angle (3 — 30 
degrees (Fig. 3), a = 45 degrees, sensor offset distance is nine pixels, and 
the distance moved per step, 6x is one pixel. There is significant interplay 
between the sensor arm angle (/?) and the rotation angle (a) parameters. 
When the two parameters are equal (for example, both 45 degrees), the effect 
on network formation is to contract the initial network that emerges. If a < P 
the contraction behaviour is increased. If, however, a > f3 then spontaneous 
branching will appear during the network formation. The branching occurs 
because the larger rotation angle places an agent's sensory apparatus outside 
the zone of the trail that the agent was following. The a/P parameters can 
be used to tune the network behaviour towards minimisation or towards high 
connectivity. After the sensory stage every agent attempts to move forward 
in its current direction (represented by an internal state from 0-360 degrees). 
If the new site is not occupied, the agent moves to the new site and deposits 
trail onto the trail map at the new location (deposition value 5 units). If the 
agent cannot move forwards successfully then no deposition is made to the 
trail map. 

Inertial movement for each agent is provided by maintaining a fioating point 
representation of the current position, as well as the discrete position corre- 
sponding to the image structure. This effectively allows the agents to 'slide 
past' one another even when the next cell is occupied (the occupation can 
only actually happen when a cell becomes free). The inertial behaviour results 
in the emergence of surging movements in the population and corresponds to 
the spatial oscillations seen in Physarum. The strength of the oscillations can 
be reduced by firing a 'change direction' event with probability 0.05 (in these 
results). When the event is triggered for an agent, the fioating point position is 
restored to the discrete position and a new direction is randomly selected. The 
effect is to dampen the surge of movement caused by the inertial behaviour. 

The list comprising the agent population is traversed in a random order for 
both the sensory and movement stages to avoid any infiuence from sequential 
positional updates. In the current version of the model the maximum popula- 
tion size is manually assigned, based on a proportion of the size of the problem 
configuration area. Work is in progress to automatically assign the population 
size in relation to food source availability. 

The experiments on the model Plasmodium were designed and configured to 
refiect the 'wet' experiments reported in this report. The results images show 
the experimental field within the confines of the 'dish'. The food sources corre- 
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(a) t=0 h 



(b) t=12 h 




(c) t=24 h 
Fig. 4. Fusing of active zones. 



spend to the oat flake positions used in the 'wet' experiments and are indicated 
as smaU dark spots. The position of the Plasmodium is indicated as the rough 
textured masses. Evolution of the model for each experiments proceeds from 
left to right. 



5 Results: Manipulating active zones 

The operation Fuse(^i, yl2), fusing of two active zones Ai and A2 can be 
implemented via collision of the active zones. An example is demonstrated in 
Fig. 4. 

Five oat flakes were arranged in a line on a humid fllter paper. Two pieces of 
Plasmodium were placed nearby the extreme flakes, one piece near South- West 
flake, another piece near North-East. The active zones are formed as pseudopo- 
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(a) t = 63 steps (b) t = 105 steps 




(c) t = 147 steps (d) t = 294 steps 



Fig. 5. Evolution of model, simulating fusing of active zones. 

dia and propagate towards the center of the chain (Fig. 4ab). When the active 
zones Ai and A2 collide they fuse and annihilate, Fuse(Ai, A2) = (Fig. 4c). 
Depending on the particular circumstances the new active zone (the result 
of fusing) may become inactive, transform to protoplasmic tubes, or remain 
active, Fuse(74i, ^2) = A, and continue propagation in a new direction. 

The model output (Fig. 5) shows the inoculation and multiplication of the 
Plasmodium at the outer nodes. The growing plasmodium is attracted to the 
inner food sources and the propagation continues inwards from each direction. 
The Plasmodium fuses and maintains its surface area spanning the array of 
nodes. The model plasmodium was seen to periodically oscillate it's position 
over the array whilst maintaining coverage. 

Multiplication of an active zone, the operation MuLT(yl) = {^1,^2}, can be 
implemented by placing sources of nutrients nearby the protoplasmic tubes, 
or inactive zones (Fig. 6). An example is shown in Fig. 6. A chain of oat flakes 
is connected by protoplasmic veins, we add two new oat flakes to evoke new 
active zones (Fig. 6a). Ten hours later two active zones Ai and A2 are formed. 
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(a) t =0 (b) t =10 h 

Fig. 6. Multiplication of active zone. 




(a) t = 63 steps (b) t = 126 steps 




(c) t = 147 steps (d) t = 273 steps 

Fig. 7. Evolution of model, simulating multiplication of active zones. 

each pseudopodium travels to its unique oat flake (Fig. 6b). 
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(a) t =0 h (b) t =12 h 



Fig. 8. Illustration of how active zone can be relocated to a domain without food 
sources, operation Direct(74, v). 

The model output initially recreates the spanning tree from the outer nodes, 
using the method in (Fig. 5). Two more food sources are added to either side 
of the array and the Plasmodium sends two streams outwards to engulf the 
sources (Fig. 7). When the food sources have been engulfed, the Plasmodium 
shifts in position by redistributing its component parts to cover the (diamond 
shaped) area created by the addition of the two new nodes. 

In some case we may need to translate active zone not to another graph node 
(source of nutrients) but to a domain of an active space. Moreover, it may be 
necessary to provide the active zone A with certain initial velocity vector v to 
the zone continues its propagation in the predetermined direction. Such oper- 
ation Direct(A, v) is executed using additional sources of nutrients (Fig. 8). 
Given a chain of oat flakes connected by protoplasmic tubes we added two new 
flakes on the west of the chain (Fig. 8a). An active pseudopodium, represent- 
ing the active zone A sprouted from the 'old' flake lying between projections 
of the new flakes. It continues propagating along the bisector separating two 
new flakes. The active zone continues its propagation towards West till collide 
with the wall of the Petri dish. Meantime, two more active zones are formed 
to connect new oat flakes to the existing protoplasmic network (Fig. 8b). 

The model was not able replicate the experimental results of Fig. 8. As soon 
as a food stimulus is removed from the environment the model Plasmodium 
begins to shrink. We speculate that the agent population is responding too 
quickly to the changing environment, and so the inertia of the agent population 
(in the current version of the model) is not as strong as observed in the real 
Plasmodium. The model of Physarum machine requires some direct stimulus 
to move the Plasmodium active zone. If it is not possible to use food source 
stimuli to move the active zone then it may be possible to use a synthetic 
hazard (such as the response of real plasmodia to certain frequencies of visible 
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(a) t =0 h (b) t =11 h 

Fig. 9. Major restructuring of the graph. 
Hght) to drive the active front to a new position by repulsion. 

Translation of active zone may lead to, or used for, undertaking of the major 
restructuring of the data storage graph. Thus, a parallel partial shift of the 
graph chain is demonstrated in Fig. 9. Three new oat flakes are places alongside 
the chain of oat flakes already connected by protoplasmic tubes (Fig. 9a). In 
few hours the Plasmodium occupies new oat flakes. Moreover, part of the 
protoplasmic graph coaligned with new oat flakes is shifted to new oat flakes 
(Fig. 9b). Abandoned protoplasmic veins, former edges of the shift part, are 
visible as white tubes. 

Such a major restructuring of the graph was caused by relocation of the active 
zone. At first the active zone travelled North-East (Fig. 9a). Addition of new 
oat fiakes caused the zone to switch to the new location and move South- West 
(Fig. 9b). 

A major restructuring was also observed in experiments with the particle 
model. Note that the particle model exhibits distinct differences to the actual 
organism. The graph restructuring with the real Plasmodium is affected by 
the remnants of old protoplasmic veins. Such debris is not present in the 
model (primarily since the 'veins' in the model merely consist of protoplasmic 
fiow itself and have no structural features). The shift in active front position 
is dramatically illustrated in Fig. 10a where the front moves from its central 
position to engulf the three new nodes on the left. Fig. 10b continues the same 
experiment, deleting the three leftmost nodes and adding three new right-side 
nodes, resulting in a significant surge of the active front to the right side of 
the 'dish'. 

In many examples above we observed formation of protoplasmic tubes, con- 
necting two geographically closest sources of nutrients. Given a chain of oat 
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(a) t=1372, 1428, 1484, 1568, 1820 steps 




(b) t=2016, 2156, 2212, 2296, 2352 steps 



Fig. 10. Major restructuring of the graph in the particle model, (a) Active front 
shifts towards new stimuli on the left, (b) Active front shifts towards new stimuli 
on the right. 

flakes, an additional oat flake is added to the experimental container. What is 
exact mechanism of inclusion the new flake in the graph? Do all nodes of the 
graph develop active zones, which travel to the flake while competition with 
each other for the flake? 

Our experimental observations show that if the graph remains connected then 
only one active zone, heading for the new oat flake, emerges. This happens 
due to sychronization of activity in the whole protoplasmic network. 

This may indicate that Plasmodium flrst 'decides' which part of the proto- 
plasmic graph closest to the recently added source of nutrients and only then 
generates an active zone. All parts of the Physarum machine sense the chemo- 
attractants coming from the new source of nutrients, however the node closest 
to the new source somehow suppresses, or inhibits, activity of other nodes. 
Existence of connections between nodes of the Physarum machine is a pre- 
requisite for the inhibition. 

To prove the point we arranged oat flakes in a chain and allowed the chain 
to be spanned by the Plasmodium (Fig. 11a). When the flakes are connected 
by protoplasmic tubes we add new oat flake, the eastmost flake on the pic- 
ture (Fig. 11a). To break communication between nodes of the Physarum 
machine we cut through protoplasmic tubes connecting three Northern nodes, 
cut places are shown by lines in Fig. 11a). Due to breakup in communication 
all nodes of the plasmodium reacts, almost simultaneously, to the addition 
of new source of nutrients by sprouting pseudopodia (Fig. lib). In few hours 
cut tubes are restored (Fig. 11c). This reinstates communication between all 
parts of the Physarum machine. The distant, from the new flake, node of 
the Physarum machine des-activate their active zones (cease propagation of 
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(c) t=17 h 

Fig. 11. Cutting edge leads up to de-synchronization of the Physarum machine. 



pesudopodia). In the result of restored communication only one active zone 
remains, and the new node is connected to the storage graph by a single edge. 



Just cutting protoplasmic tubes does not lead to formation of new pseudopodia 
(Fig. 12). When tubes are cut they do usually fuse back again in few hours 
without any adverse effect to the Plasmodium behaviour. 



The model Plasmodium was not able to replicate the experiment of Fig. 11, 
once again suggesting that the model is currently lacking in the mechanism 
which couples the Plasmodium behaviour across food sources. 
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(a) t=0 h (b) t=14 h 

Fig. 12. Absence of adverse reaction to tube cutting. 
6 Summary 



We studied the Physarum machine, an implementation of general-purpose 
Kolmogorov-Uspensky machine (KUM) in vegetative stage, Plasmodium, of 
Physarum polycephalum. Our present paper is mainly concerned with manip- 
ulating active zones (actively growing pseudopodia of the Plasmodium) be- 
cause they are main computational units of the storage modification machine, 
implemented in the Plasmodium. In laboratory experiments and computer 
simulation of the Plasmodium we executed basic operations with active zones. 
We shown how merge two active zones, to multiply an active zone, to translate 
active zone to a new node of the storage structure, to direct active zone in the 
space not occupied by data nodes. We envisage our experimental and theoret- 
ical findings will be employed in future programming of spatially distributed 
biological computers. 
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